Quantum pattern recognition with liquid-state nuclear magnetic resonance 
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A novel quantum pattern recognition scheme is presented, which combines the idea of a classic 
Hopfield neural network with adiabatic quantum computation. Both the input and the memorized 
patterns are represented by means of the problem Hamiltonian. In contrast to classic neural net- 
works, the algorithm can return a quantum superposition of multiple recognized patterns. A proof of 
principle for the algorithm for two qubits is provided using a liquid state NMR quantum computer. 
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I. INTRODUCTION 

The framework of Natural Computing tries to model 
architectures found in Nature and to apply them to a di- 
versity of computational tasks. Its biological domain is 
represented among others by neural networks and evo- 
lutionary computation. In physics, novel algorithmic 
schemes have been investigated in the context of quan- 
tum computation. Based on these ideas, we investigate 
how pattern recognition can be implemented using prin- 
ciples of quantum computing. 

Artificial neural networks involve interacting units 
called neurons designed according to neural structures of 
a brain which cooperate in order to process information 
[l[ . Similar to the brain, artificial neural networks are ca- 
pable of performing cognitive tasks such as pattern recog- 
nition and associative memory. Pattern recognition pro- 
cesses input data usually on the basis of a priori knowl- 
edge in the form of a set of memorized patterns. An input 
pattern is then classified to whichever one of the memory 
patterns it most closely resembles. Typical applications 
for pattern recognition include automatic recognition of 
objects and patterns in digital image analysis as well as 
voice and speech recognition. An associative memory is 
a system that completes partially known input pattern 
based on the stored content. Typical applications for as- 
sociative memory include content-addressable memory as 
a special type of computer memory and database engines. 

A classic theoretical model to perform the described 
tasks is a Hopfield network Q , a recurrent neural network 
with symmetric connections between individual neurons. 
Full connectivity is provided with every neuron i being 
able to interact with any other neuron j 7^ « by means 
of a response function Vi (t) = J2j^i '^ij (*) aggregated 
as a weighted sum over the current states of all other 
bipolar neurons Sj{t) = ±1 which correspond to biolog- 
ical states of not firing and firing electrical signals to its 
neighbor neurons. An activation function fi{x) of the 
neuron i then evaluates in order to define if the neu- 
ron remains in its currents state or a state flip is applied, 
e.g., Si{t + At) = ft{ri{t)), where fi{x) = sgn{x) is the 



sign function. 

For a set of p patterns V ^ {S,^ , . . . = {^'"j with 
fi ~ 1, . . . ,p and bipolar ^'^ — ±1 different definitions for 
the synaptic connection strength Wij between neurons i 
and j are possible, e.g. by the (discrete) Hebbian matrix 



1 

N 



p 



(1) 



which stores the memory information in a distributed 
manner. The input binary vector ^'"^ is imposed on the 
Hopfield network as its initial state. The dynamics of the 
system is designed such as to minimize a cost function 



(2) 



by converting the state of the network to a stable config- 
uration which represents the recognized pattern. Thus, 
the physical background of the Hopfield model offers an 
approach for quantization of the neural network. 

Quantum computers utilize special characteristics of 
quantum systems such as superposition and entangle- 
ment 3. The corresponding algorithms provide compu- 
tational effectiveness in comparison with classic routines 
for problems such as search [J] and factorization 5] . Dif- 
ferent hardware designs have been proposed for build- 
ing a quantum computer, including liquid state NMR 
[1, 0, H [§| and ion traps [13] ■ The extensive technolog- 
ical development of NMR spectroscopy in the last five 
decades made it possible for NMR quantum computing 
to become a suitable test ground for novel quantum al- 
gorithms. 

A systematization of quantum computing using adia- 
batic evolution leads to the concept of Adiabatic Quan- 
tum Computation (AQC) [llj that has been proven to 
be equivalent to the standard network model of quantum 
computation [l^ . Compared to the abstract computa- 
tional language of the quantum networks model, AQC 
often provides a more direct translation to experimental 
quantum computing. 
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The protocol for AQC is represented by a controlled 
Hamiltonian path 7i(s) with s = y G [0, 1] and run- 
ning time T . The computational problem is encoded in 
the final Hamiltonian 7i(l) — Tip and the solution of 
the computational problem corresponds to finding the 
ground state of Tip. On the other side, the initial state 
of the quantum system |'!/'(0)) is chosen to be an easily 
preparable ground state of 7i(0) — Tii. If the Hamilto- 
nian is driven from Tii to Tip slowly enough, i.e., if the 
total calculation time T is chosen to be long enough, then 
in the adiabatic limit the final state of the quantum sys- 
tem comes arbitrarily close to the problem ground state. 
AQC can also be formulated in the framework of Quan- 
tum Annealing [l^] utilizing a Hamiltonian trajectory 

Ti{s)^ K{s)Ti^+Tip (3) 

with A(0) large enough to make Tii the dominating 
term and A(l) = which reduces the Hamiltonian to 
the constantly present contribution by Tip. The trans- 
verse Hamiltonian Tii represents the driver for the adi- 
abatic evolution whereas the initial state |i/'(0)) of the 
annealing process is chosen so it corresponds approxi- 
mately to an energy eigenstate of the initial Hamiltonian 

H(0) = AjnaxTii + Tip. 

In classic terms, AQC represents minimization of an 
energy cost function. This offers a link from (usually 
irreversible) classic neural dynamics to quantum dynam- 
ics governed by unitary, reversible evolution. However, 
the classic Hopfleld network approach performs pattern 
recognition by local optimization. In the following, we 
utilize AQC which performs global optimization in form 
of ground state approximation. 

II. THEORY 

We consider a quantum neural network with N neurons 
consisting of bipolar neural states —1 and -1-1 to be repre- 
sented by a quantum system with N qubits in states |0) 
and Whereas a classic artificial neuron can assume 
only a single "fire" or "not fire" state at once, quantum 
neurons allow superpositions of these two states in the 
form of a|0) + /3|1) with \a\^ + \f3\^ = 1. 

In order to apply AQC to the computational task of 
(associative) pattern recognition, we encode the complete 
problem in the Hamiltonian 

Tip = Timem + TTirnp (4) 

with Timem representing the knowledge about stored pat- 
terns, Tiinp representing the computational input and an 
appropriate weight factor F > (see Fig. [1]). 

Whereas the classic Hopfield network utilizes dynam- 
ics defined by the memory patterns alone, the dynamics 
according to Eq. (|4]) depends on both memory and input 
information (see Fig. [J). 

Furthermore, the encoding of the final adiabatic 
Hamiltonian as a sum of two independent terms offers 
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FIG. 1: Illustrative representation of the energy landscapes 
for the memory Hamiltonian Timem (A), the weighted input 
Hamiltonian VTLinp (B) and the total problem Hamiltonian 
Tip (C). Dashed circles represent memory patterns stored as 
energy minima of Timem', crossed dashed circle represents the 
input pattern as energy minimum of Tiinp- 



an advantage with respect to the practical implementa- 
tion of the process. Using appropriate approximations 
(such as Suzuki- Trotter), the computational evolution of 
the system can be approximated by separate periods of 
unitary evolution each according to the individual terms 
building the total instantaneous Hamiltonian. This im- 
plies that the memory precompilation, i.e., the trans- 
lation of dynamics according to Timem, into elementary 
physical interactions such as pulse sequences in NMR 
quantum computing systems has to be executed only 
once for an individual memory set. The generated low- 
level commands can then be reused for variable input in- 
formation. The dynamics according to the input Hamil- 
tonian Tiinp can be obtained and imposed independently. 

In direct analogy to energy function ([2]), the memory 
Hamiltonian can be defined by the coupling strengths 
between the qubits 

Timem = ^ 2 X] '^ij'^i'^h 

where af is the Pauli z matrix on qubit i and Wij is 
the Hebbian matrix ([1]). A more general choice of the 
memory Hamiltonian is discussed in Appendix \^ 

For the retrieval Hamiltonian Tiinp we consider a single 
input pattern of length N . For associative memory 
applications we can extend a noncomplete input vector 
^inp q£ igQg^}j n < TV by setting the values for N — n 
unknown states to zero. In contrast to the classic neural 
network, we impose the input pattern on the dy- 
namics of the system by an additional Hamiltonian term 
Tiinp which can be defined as 

w™p = Erv- (6) 

i 
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FIG. 2: Illustrative representation of classic and quantum dy- 
namics within the energy landscape. Classic dynamics drives 
the initial state (crossed circle) which represents the input 
pattern toward the final state (solid circle) which represents 
a nearby memory pattern. Quantum dynamics transforms 
the energy landscape from the initial Hamiltonian close to 
Hi (A) via intermediate Hamiltonian Ti.{t) (B) towards the 
final Hamiltonian equivalent to the problem Hamiltonian Tip 
(C). The current state of the system remains the ground state 
of the instantaneous Hamiltonian with the "empty memory" 
state \tpun) (0 for t = (represented by the star in A) and 
the final state of the system corresponding to the outcome 
memory pattern for t = T (represented by the solid circle in 
C). 

For a given state |^) representing a pattern ^, the energy 
corresponding to Hinp is given by £'™P(f) = —n + 2h and 
depends on the Hamming distance h — 0, . . . ,n between 
^inp and ^ = [^1, . . . , ^„], i.e., the number of positions for 
which the corresponding entries are different. The ex- 
ternal field defined by Timp thus creates a scalar metric 
proportional to the Hamming distance between and 
memory patterns. This permits a quantitative compari- 
son between input and memory patterns by shifting the 
energy levels of the memory Hamiltonian TLmem so the 
more similar memory patterns have lower energy com- 
pared to the alternative patterns. It can be shown that 
for the case of a single stored pattern [p = 1), the up- 
per bound for the weight factor in Eq. ([4]) is given by 
r < 1 — ^ which ensures that the correct memory pat- 
tern is the ground state of the system (cf . Appendix [B| . 

Finally, we define the initial conditions of the AQC pro- 
tocol. In comparison with the classic Hopfield network, 
the initial state of quantum pattern recognition systems 
does not correspond to the input pattern but is chosen 
as 

1^(0)) ^ \^Pun) = -i E 1^) 

representing an "empty memory" with uniformly dis- 
tributed probability for all possible state configurations. 



As corresponding initial Hamiltonian, we choose an eas- 
ily constructible sum over single qubit operators Tii = 
i^-(l — erf) which can be implemented experi- 
mentally in a straightforward way. Since the problem 
Hamiltonian Tip ([4]) is diagonal, the transverse form of 
TLi ensures the transformation of the initial state \^un) 
([T]) toward the outcome memory pattern. 

III. EXPERIMENTS 

For the simulations and experiments, we considered a 
heteronuclear two spin system corresponding to the ^H 
and ^^C nuclear spins ^ of ^'^C labeled sodium formate 
dissolved in D2O. The experiments were performed at a 
temperature of 27°C using a Bruker AC 200 spectrom- 
eter operating at a Larmor frequency of 200 MHz for 
iH and 50 MHz for "C. The heteronuclear ^H-i^C cou- 
pling constant is J = 195 Hz and the relaxation times 
are Ti{^YL) = 1.6 s, Ti{^^C) = 2.7 s, T2{^Y^) = 130 ms, 
and T2(^'^C) = 60 ms. In this section, we use normalized 
Pauli operators as conventional in NMR [isj . 

In order to provide an experimental demonstration for 
the concept of quantum pattern recognition, we imple- 
ment a basic two-neurons Hopfield network using a liq- 
uid state NMR quantum computer. The single connec- 
tion w = W12 — W21 between two qubits corresponds for 
w = — 1 to patterns [— 1,+1] and [-1-1,-1] with differ- 
ent bit values whereas w = +1 corresponds to patterns 
[—1, —1] and [-1-1, -1-1] with identical bit values. Following 
Eq. ([4]), the problem Hamiltonian can be written as 

Tip = -walal + T [^''^t + ^^'^2) ■ (8) 
The NMR free evolution Hamiltonian 

Tip = 2t: Juh^c + '^t^vhcJh + 2t:vc(Jc (9) 

for the considered molecule shows a constant positive 
coupling. The variable memory weight w according to 
Eq. ([5]) can be implemented by changing the sign of the 
effective coupling Hamiltonian using additional pulses 
[isL IT9I ]. In the case of two qubits considered here an 
alternative approach which avoids the need for addi- 
tional pulses is the following. We can map the problem 
by defining the offset frequencies vh = — wFJ^™^ and 
vc — —wTJ^™^ and redefining the initial Hamiltonian 
as problem dependent TL[ = 2TTw{afj + a^). The ini- 
tial quantum state li^un) of Eq. ([7]) corresponds to the 
ground state of TCi for w = —1 and to its highest excited 
energy eigenstatc for w = +1. Our choice of vh and vc 
consequently encodes the desired \ipres{w,£,^"P)) as the 
ground state, respectively, the corresponding highest ex- 
cited state of the Hamiltonian TL{T). In classic terms, 
the first case represents minimization of a cost function 
whereas the second involves its maximization. 

Following Eq. ([7]), the desired initial state for the 
considered heteronuclear two spin system is \^un) = 



5(|00) + |01) + |10) + |11)) with corresponding density 
matrix p„„ = \'ipun){''Pun\- Rewriting this density ma- 
trix in terms of spin operators, ignoring the identity 
term and proportionahty constants, the initial state can 
be represented by the simphfied and traceless operator 
p(0) = a'fj+aQ + 2afj(TQ. This state can be created start- 
ing from the thermal equilibrium using standard methods 
based on spatial averaging 20]. 

The transverse form of the initial Hamiltonian Tii sug- 
gests that pulse sequences applied on the transverse di- 
rection X can be used not only to generate but also to 
effectively control the initial Hamiltonian. More specif- 
ically, Tii is represented by a transverse radio frequency 
(rf) field in positive or negative (depending on w) x di- 
rection on both spins. The amplitude of the rf field A(t) 
is reduced linearly from A(0) = A^ax to A(r) = during 
the fixed evolution time T = 50ms which is shorter than 
relaxation times Ti and T2. The offsets vh and vc in 
Eq. ^ were set to values ±100 Hz and Hz correspond- 
ing to values = ±1 and S^i — Q, respectively. 

In the simulations and experiments of quantum anneal- 
ing, the time evolution of the system was discretized into 
L steps of length At — T/L with constant Hamiltonian 
7i(i) for the time interval [<, At]. The continuous evo- 
lution is approximated for L — )■ 00 and At — > 0. Our nu- 
merical results suggested that L = 100 leads to sufficient 
accuracy. We chose F = 0.5 for the weight parameter of 
the retrieval Hamiltonian term. We set Amax = 600 Hz, 
where in simulations we get a reasonable value of 98% for 
the overlap tr (pfPexp) between final state of the system 
Pf and expected outcome state 

Pexp • 

The result of the quantum annealing is read out by 
applying ^ pulses followed by detection of the ^H or ^^C 
signals. The ^H signals were acquired as single scan while 
^^C signals were acquired using 16 scans. The spectra 
recorded after the read out are shown in Fig. [31 A rea- 
sonable match between the experimental results and sim- 
ulations (data not shown) is found. The first and third 
columns represent ^H spectra which encode the computa- 
tional results as presented in Table |T1 The first four rows 
represent applications of associative pattern recognition 
with n = 1 and N = 2 where the algorithm evaluates 
the Hamming distance between the input and the partial 
memory patterns. A single peak at frequency vh^^ cor- 
responds to a unique recognized pattern. The peak's po- 
sition and sign direction identify the pattern recognition 
results, e.g., left negative peak identifies the state | — 1, 1) 
corresponding to the output pattern [-1,1]. For the case 
of the blank input pattern = ^2"'° = Oi the quantum 
neural network returns its memory content as consistent 
patterns [-1,1], [1,-1] for w — ~1 and [-1,-1], [1,1] for 
w = +1. The resulting quantum state represents a su- 
perposition of pattern states identified by corresponding 
spectral peaks. For completeness, we show in the second 
and fourth column of Fig. [3] corresponding ^^C spectra. 
Peaks a.t Vc ^ ^ provide redundant information for the 
identification of the patterns. 
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FIG. 3: Experimental and ^^C spectra for five combina- 
tions of ^l"^ and ^2"'' two possible values of w, where 
w = —1 and w = +1 correspond to stored patterns with op- 
posite and equal bit values, respectively. The spectral range 
in each spectrum is ±240 Hz. 
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TABLE L The results of the pattern recognition represented 
by the corresponding pure quantum states. 



IV. CONCLUSION 

In summary, we have developed a new theoretical ap- 
proach for quantum pattern recognition and implemented 
it using liquid state NMR techniques in a simple exam- 
ple with two qubits. In the case of an associative mem- 
ory, an incomplete input typically yields several memory 
patterns which are equally similar to this input. In con- 
trast to classic neural networks, a quantum neural regis- 
ter can represent a superposition of recognized patterns. 
Whereas linearly combined classic mixture states, e.g., 
of the form sgn(±^^ ± ± ^'") 1] do not provide direct 
information about corresponding memory patterns, the 
quantum superposition states allow unambiguous read- 
out of distinct patterns which form the superposition by 
means of an appropriate (iterative) measurement proce- 
dure. 

The projector memory Hamiltonian enables the mea- 
surement of similarity between the input pattern and 
stored patterns as well as allows a differentiation between 
memory patterns and their reversed patterns — in con- 
trast to classic Hopfield network (cf. Appendix Since 
the problem is directly encoded into the final Hamilto- 
nian Tip, no ancilla qubits are required. Furthermore, 
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AQC dynamics rules out closed state cycles of the form 
S{t + T) — S(t) with a period T > 1 which are possible 
for classical neural networks [H. 

In the experiments, we realized the desired morph- 
ing between the Hamiltonians using quantum annealing. 
This allowed us to exploit the experimentally available 
control terms and did not require an iterative synthesis 
of the effective Hamiltonian trajectory [2^ . Furthermore, 
we performed both ground and excited state AQC. 

Several characteristics of the introduced concept re- 
main open and require further research. The compari- 
son of maximum storage capacity of the presented AQC 
based model and classic Hopfield network is an open 
problem that motivates additional study. Furthermore, 
time requirements for pattern retrieval can be investi- 
gated by studying the minimum energy gap during the 
course of the evolution [23| as well as employing heuristics 
related to success probability Future investigations 
could also involve optimized annealing procedures which 
could provide significant improvements over alternative 
schemes [l^. An interesting enhancement of the pro- 
posed scheme could be introduced by considering asym- 



metric synaptic weights in neural modeling [17[. Future 
directions of research also include the extensions to larger 
spin systems using liquid state NMR quantum computa- 
tion with the control of individual couplings of the ef- 
fective Hamiltonian Ti.{t) executed by means of pulse se- 
quences and refocusing schemes [3, [U, [2^, [2a|- The 
developed techniques can be potentially applied to alter- 
native hardware designs, such as ion traps with magnetic 
field gradients [l^l which have been already contextual- 
ized in quantum neural research [2^ and superconducting 
circuits [29j . 
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APPENDIX A: PROJECTOR MEMORY 
HAMILTONIAN 

In addition to the definition of the memory Hamilto- 
nian of Eq. ([S]) , we introduce an alternative encoding of 
memory information based on a generalization of quan- 
tum search projection operators [2l[. It can be defined 
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as a memory Hamiltonian 



(Al) 



which is diagonal in the computational basis. The defini- 
tion represents the ideal case of memorized patterns en- 
coded in the (degenerate) ground state energy of Ti-mem 
whereas not-memorized patterns correspond to excited 
states. In combination with the input Hamiltonian Timp 
([6]), the bound for the input weight is given by F < ^ 
for arbitrary values of p (cf . Appendix [B]) . 

Alternatively, the memory Hamiltonian can be speci- 
fied as a projector on the memory state 



'^mejn ^ Is 

with memory state 



niem\ / /-jnem 



(A2) 



(A3) 



In the case of both projector memory Hamiltonians, the 
initial Hamiltonian can be chosen as a projector 7i,; — 

un I ■ 

Using the definition (|A2p . we can solve the prob- 

lem of quantifying the level of similarity between the in- 
put pattern and all memory patterns. For this, the 
probability for measuring the memory pattern is inter- 
preted as a measurement of relevance. For comparison, 
the initial probabilities for the memory state (|A3|1 are 
equally distributed with = ^ 1, com- 

bination with the input Hamiltonian Timp ([6]) , the ground 
state of the problem Hamiltonian Tip = 'H^mem + ^"Hmp 
is 



1 



1 - 2rSh, 



m + o(r2), (A4) 



where = d^iC"^ jS,^) is the Hamming distance be- 
tween i^** — [Ci , ■ • ■ ,S,n] ^iid the retrieval pattern, Sh^^ = 
hf_i — (h) is the deviation of the Hamming distance from 
the average relative Hamming distance (h) = i X]^=i ^t^ - 
Thus, the input Hamiltonian shifts equally distributed 
weights in the memory state (|A3[) according to the cor- 
responding Hamming distances. 

We outline this effect using perturbation theory ^3d\ for 
the problem Hamiltonian Ti.^ in relation to the parame- 
ter r. The ground state energy and the corresponding 
eigenstate can be represented by 



(1) 



F'lK 



0(f3), 



(2) 



0(f3).(A5) 



The first order approximation energy is given by E'q — 
(^^rnem^j^znp^-^rnemy ysiug the definition of the input 

Hamiltonian, it follows for = [<^f , . . . , £,„]■ 



(A6) 



where s^, = ^'"^-^'^ is the scalar product between vectors. 
Since = n — 2h^ with the Hamming distance = 
c?h(C™^,?'^) for the first n elements, it follows 

VP ^ 

Using (^^li^") = S^f^, we obtain for the first order energy 
correction 



= _n(e"^'"ic^'") + -E(ri/iMi^" 



= -n + 2(h) 

with average relative Hamming distance (h) = ^ /i^. 
The first order approximation state is given by 



4^') - - E 



By definition of the Hamiltonian Tiji^em (|A2p , the energy 



levels are EJ^'''^ = and E]^^l[^ 



1. Since \E]^^"^) form 



an orthonormal basis, it follows using Eq. (jA6 



\E 



(i-ir"")(r'="i)^E^^'i^' 



with average relative scalar product (s) = i = 
n — 2{h). Using state approximation (lASp . we obtain for 



the zero and first order 



\EE) = \E, 
1 

VP 



/ 



F|^, 



We rewrite this expression in terms of the Hamming dis- 
tance, thus obtaining expression (|A4[) . 



1 - 2T5h, 



1^^) + o(f2 



In case of all memory patterns being equidistant from 
the input pattern with regard to the first n positions, the 
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product basis 

FIG. 4: Final measurement probabilities according to mem- 
ory set (|A7|) and input pattern (|A8|) . The peaks corre- 
spond to memory patterns = | — 1, —1, —1, —1, —1) = |1) 
(peak A), .^^ = | - 1,-1,-1,-1-1,-1) = |3) (peak B) and 
= I - 1, -1, +1, -1, -f 1) = |6> (peak C) in binary and dec- 
imal representation, respectively. 



deviation Sh^ equals zero so that corresponding measure- 
ment amplitudes are equally distributed. Otherwise, the 
amplitudes are shifted according to whether the distance 
to the input pattern is above- or below-average. 

For completeness, we derive the second order approxi- 
mation energy given by 



E, 



(2) _ 



E 

fc=i 



^mem 



rpmem 
Eq 



It follows E^^^ = -4A/i2 where A/i^ = {h'^) _ (h)"^ is the 
standard deviation for hi and (h^) = ^ J^i • 

In order to affirm the bias beyond the approximation to 
first order, we demonstrate the corresponding probabil- 
ity distribution for an exemplary memory set and input 
pattern 



V = 




(A7) 
(A8) 



with Hamming distances between the input pattern and 
the memory set given by c?i = dniC"^ j^^) = 0, (i2 = 1, 
and ds — 2. Fig. [K\ shows for F = O.I the distribution of 
measurement probabilities with pi « 0.476, p2 ~ 0.308, 
and P3 « 0.216. 



APPENDIX B: INPUT WEIGHT 

For p = 1, we derive the bounds for the weight factor F 
in combination with the coupling memory Hamiltonian 



'Hmem ©■ Sincc the Hamming distance between two 
patterns is symmetric under flips of the same positions, 
the memory pattern of length N and input pattern of 
length n < N can be written in the form 



N 



1 -1 



N- 



G 



1-1 



1+1 



where m is the Hamming distance between the patterns 
for the first n elements. An arbitrary vector ^ of the 
length N can thus be written in the form 



1+1 



1-1, 



1+1... +1 



{n—m) — mi 




where mi, m2, and represent the distances betweens 
f and three respective segments of the memory pattern 
of the length n—m, m, and N—n. Using the definition 
of memory and input Hamiltonians, we obtain 



E^e) 



2m), 



N 1 

+ - - F n 

2 2 ^ 



+ F (n — 2m) . 



£;""^"(-^i) + Fi;"P(-Ci) 
N 1 

- -y + 2 

EP{£_) = E""'"'{^) +TE"'P{^) 

I 
2 

F (n — 2 (mi + (m — m2))) , 



-^[iV-2Mf 



where M = mi + m2 + ma if d{£^,£^^) < d(^, —S,^) and 
M ^ N — (mi + m2 + ms), otherwise. 

Since EP{-^^)-EP{^^) cx F (n - 2m), for all valid val- 
ues of F the memory pattern closer to the input pattern 
also has lower energy level. In the following, we assume 
^mp closer to the original pattern with < m < f 
where f = [^J . For the system to return a valid an- 
swer, must be the global ground state of the problem 
Hamiltonian Hp ^ with 



EP{e)<Epi^^e)- 

From condition (|Bip . it follows 

iV^ - [TV - 2 (mi + m2 + mg)]^ 



F (m2 



mi) < 



AN 



(Bl) 



(B2) 



In order to obtain the upper bound for the value of F, we 
attempt to minimize the energy value corresponding to 



8 



^. First, we assume mi = so that the overlap between 
first segments of resp. ^'"'^ and the corresponding 
spins of ^ is maximal. Furthermore, we assume ma = 
so that the overlap between the third segment of and 
^ is maximal as well. 

Since m2 < m and m is bounded by ^ , we obtain from 
Eq. (IB2|) the condition 



r < 1 



n 

2N' 



A similar bound can be defined in combination with 
the projector memory Hamiltonian 7i"„e„j (|A1[) for an 
arbitrary value of p. We consider a subset of the memory 
set Vmin = {f^, ■ • • j^^"'"} C V, the patterns of which 



show the minimal distance to the input pattern with 
e r^,n ■■ = hrmn- The Condition for the 

weight factor F can thus be given with 



-F (n - 2/i™„) < EP{^ ^ e) = 1-r 



Since we are interested in a upper bound for F, we assume 
h = 0, i.e., ^ completes the partial input vector. With 
< hmin ^ we obtain the corresponding condition 



F < 



1 

2^' 



